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ABSTRACT 

We present a composite spectrum of 60 long duration gamma-ray burst (GRB) afterglows with 
rcdshifts in the range 0.35 < z < 6.7 observed with low resolution optical spectra. The composite 
spectrum covers the wavelength range 700-6600 A in the rest frame and has a mean signal-to-noise 
ratio of 150 per 1 A pixel and reaches a maximum of ~300 in the range 2500-3500 A. Equivalent 
widths are measured from metal absorption lines from the Lya line to ^5200 A, and associated metal 
and hydrogen lines are identified between the Lyman break and Lya line. The average transmission 
within the Lyman forest is consistent with that found along quasar lines of sight. We find a temporal 
variation in fine-structure lines when dividing the sample into bursts observed within 2 hours from 
their trigger and those observed later. Other lines in the predominantly neutral gas show variations 
too, but this is most likely a random effect caused by weighting of individual strong absorption lines 
and which mimics a temporal variation. Bursts characterized with high or low prompt GRB energy 
release produce afterglows with similar absorption line strengths, and likewise for bursts with bright 
or faint optical afterglows. Bursts defined as dark from their optical to X-ray spectral index have 
stronger absorption lines relative to the optically bright bursts. The composite spectrum has strong 
Call and Mgn absorption lines as commonly found in dusty galaxies, however, we find no evidence 
for dust or a significant molecular content based on the non-detection of diffuse interstellar bands. 
Compared to starburst galaxy spectra, the GRB composite has much stronger fine-structure lines, 
while metal absorption lines are weaker. 

Subject headings: Galaxies: high-redshift gamma rays: bursts 



1. INTRODUCTION 

Optical spectra of long duration gamma-ray burst 
(GRB) afterglows provide a means to study the physical 
condition of the interstellar medium in their host galax- 
ies through the associated metal absorption lines. While 
the host galaxies themselves are generally too faint to 
allow for detailed spectroscopic studies, spectra of the 
afterglows obtained shortly after the burst can achieve 
much higher signals. High resolution spectra of individ- 
ual bursts have shown that GRBs reside in galaxies with 
large gas column densities and sub-sol ar metallicities 
(jVreeswiik et all 12004 IChen et al.ll2005[ ). Host galaxy 
abundance ratios and metallicities have bee n measured 
(|Savaglio et al.1 120031 : IVreeswiik et~aTl 120061 ). and evi- 
dence for galaxy outflows has been established from the 
presence of mu ltiple absorption com ponents separated in 
velocity space (jSchaefer et al.ll2003l ) and from correlation 
of absorption line species with host galaxy luminosity 
(|Chen et al.ll2009l) . The distances from the burst site to 
the absorbing material has been determined to be larger 
than a few tens of pc dDessauges-Zavadsky et al.l 120061: 
iProchaska et~aTI 120061: IChen et al.1 l2007t ID'Elia et ail 
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2010), while i n other cases the time variations of fine- 
structure lines (IVreeswiik et al.l2007UD'Elia et al.ll2009at 
iLedoux et al.ll2009D or the presence of vibrationally ex- 
ited H 2 molecules (jSheffer et al.l l2009h indicated a dis- 
tance of 0.4-2 kpc. While this is in itself interesting for 
the study of the nature and physical conditions present in 
high redshift star forming galaxies, one of the goals of af- 
terglow absorption line studies is to constrain the nature 
of the progenitor of the GRB, but the diverse properties 
of GRB afterglows have made this difficult. 

At the average bur st redshift of (z) ~ 2.2 
(| Jakobsson et al.l l2006bl see also Web site 6 with an 
updated sample) o bserved with the Swift satellite 
(jGehrels et al.l [2004f ) . the rest-frame UV spectra reveal 
that GRBs occur in galaxies which contain large amounts 
of neutral gas. Generally, the hydrogen column den- 
sity determined from the GRB spectra lie in the damped 
Lyman-a (DLA) regime log (JV(Hi)/cm -2 ) > 20.3, and 
both hydrogen and metal column densities are typ- 
ically larger tha n observed in inter v ening PL As to- 
wards quasars (jSavaglio et al.l l2003t Uakobsson et al.l 
l2006al ). Quasar sight lines with DLAs intersect and 
probe the interstellar environment of a galaxy at an 
impac t parameter which is ty pically less than 10 kpc 
(e.g., iMoller et all I2002L 12004ft . although some candi- 
dates DLA galaxies are found at even larger impact 
parameters (jChristensen et all 12007) in the outskirts of 
the galaxies. In contrast, long durati on GRBs explode 
in yo ung star forming galaxies (e.g., IChristensen et al.1 
2004), and are located close to the most luminous re- 
gion in the host dBloom et al.l l2002l: iFruchter et "aT]|2006l: 
ISvensson et al.ll2010l ). The difference between the two 

6 http: //www.raunvis .hi . is/~pja/GRBsample . html 
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samples can be described by the diversity in impact 
parameters (iProchaska et al.1 [2007bl : IFvnbo et all 120081 : 
iPontzen et al.ll2010D . 

While high spectral resolution studies can reveal the 
impact on the environment near the GRB event, low 
spectral resolution studies obtained for fainter afterglows 
also provide valuable information. However, significant 
discrepancies have been found between GRB afterglows 
where high resolution spectra seem to indicate lower 
metallic ities than the aftergl ows observed with lower res- 
olution ()Ledoux et al.l 12009), although the effects or bi- 
ases which cause this discrepancy are not understood. 

The current knowledge of the interstellar medium 
(ISM) in GRB host galaxies relies on a few well studied 
afterglows at a range of redshifts. Additional knowledge 
can be gained by taking advantage of a large set of GRB 
afterglows at different redshifts which sample different 
rest-frame wavelengths. By creating a composite spec- 
trum, we can study the typical properties in a high qual- 
ity spectrum over a large rest-frame wavelength range. 
Here we present a composite spectrum of 60 GRB after- 
glows from the largest s ample which has been an alyzed 
in a consistent manner (TFvnbo et al. l l20rH 1201(1 . The 
composite spectrum allows us to measure which lines are 
typically present in a GRB afterglow and their average 
equivalent widths and metallicities. 

This paper is organized as follows. In Sect. [5] we ex- 
plain the benefits and limitations for the physical inter- 
pretation of the composite afterglow spectrum. We de- 
scribe the full sample and how the composite spectrum 
is created in Sect. [3] In Sect, fj] the global properties are 
investigated and the list of absorption lines is presented, 
while in Sect. [3] we investigate specific absorption lines 
in detail. In Sect. [5] the absorption line strengths are 
compared with other high redshift galaxies. Metallicities 
are derived in Sect. [3 We discuss the results and po- 
tential bias effects in the interpretation of the composite 
spectrum as well as future prospects in Sect. [5J 

2. MOTIVATION FOR CONSTRUCTING A COMPOSITE 
SPECTRUM 

The goal of this study is to use a unique database of 
low-resolution spectra to explore whether fainter lines are 
detectable when combining all spectra compared to de- 
tections in individual, lower signal-to- noise (S/N) spec- 
tra. Since all afterglow spectra seem to be unique, we also 
want to explore specific absorption line features and their 
possible dependence on global properties for the GRB ex- 
plosion. A detailed investigation of the variations of each 
detected absorption line in the set of spectra will be pre- 
sented in a forthcoming paper (A. de Ugarte Postigo et 
al., 2011, in preparation). 

A potential benefit of having a high signal-to-noise ra- 
tio afterglow composite spectrum which covers a large 
rest-frame wavelength range is to compare existing or 
future individual low-resolution afterglow spectra to an- 
swer a set of questions: (i) Does the (new) optical af- 
terglow spectrum have unusual absorption lines or ratios 
thereof? (ii) What absorption lines are detectable given 
the GRB redshift and the S/N of the data? (iii) Is it pos- 
sible to derive redshift of a very noisy afterglow spectrum 
by convolving it with the composite spectrum? 

Certain limitations exist for the analysis of a composite 
spectrum too, where some quantities are better explored 



with high resolution optical spectra. First, the velocity 
spreads of individual absorption lines are not possible to 
analyze at low spectral resolution. Second, some lines 
may be particularly blended with other lines, for exam- 
ple molecular H2 lines which fall in the Lya forest. In 
addition, iron fine-structure lines are affected by blending 
in the low resolution data, so a detailed analysis of their 
temporal variations has to use higher spectral resolution 
data. 

Individual afterglows have shown a wide range of 
metallicities in GRB hosts, and the composite is biased 
towards the brightest afterglows in the full sample as ex- 
plained in Sect. [3] However, with this in mind we can 
use the wide rest-frame wavelength coverage of the com- 
posite spectrum combined with the detection of a large 
number of weak absorption lines, to derive the average 
column density of a large number of metal species in 
Sect. [71 This presents a unique possibility for the study 
of the interstellar medium (ISM) in faint, high redshift 
galaxies, although the diversity of metallicities probed 
by GRBs have to be investigated using individual higher 
resolution spectra. 

3. AFTERGLOW SAMPLE 

The largest systematic analysis of low resolution opti- 
cal spectra of long durati on GRB afte r glows performed 
to date is presented in IFvnbo et al.1 (|20 09. hereafter 
F09). While the majority of spectra in F09 involves 
low-resolution spectra, they also include five high resolu- 
tion spectra observed either with UVES on the VLT or 
HIRES at the Keck Telescope. To date, optical spectra 
of GRB afterglows are mostly obtained at low resolution 
because afterglows are typically faint (R > 20 mag) at 
the time of the start of the observations. The low spec- 
tral resolution data from the sample in F09 provide the 
basis for this present follow-up study. 

3.1. Sample selection 

The aim of the sample selection is to construct a sample 
of long GRBs that is selected independent of the optical 
properties of the afterglows and at the same time has 
as high completeness in high quality optical follow-up as 
possible. 

To get a well selected sample which is independent of 
the optical properties of the afterglow, and to achieve 
a high degree of completeness and good quality optical 
spectra, a set of selection cri teria was used accordin g 
to the methods described in Uakobsson et al.l ((2006b): 
(i) the burst duration is longer than 2 sec, (ii) accu- 
rate localization obtained with X-ray observations, (iii) 
small foreground Galactic extinction (Ay < 0.5 mag) 
in the direction of the burst, (iv) declination useful 
(—70° < S < +70°) for optical follow up observations, 
and finally (v) a solar distance greater than 55° at the 
time of the GRB. A total number of 146 Swift bursts 
observed between March 2005 and September 2008 fulfill 
these criteria, and optical spectra have been obtained for 
69 of these (see F09 for details on each individual object). 

From the full sample of spectra, we here focus on the 66 
afterglows which were observed with low-resolution opti- 
cal spectroscopy using the following instruments: FORS1 
and FORS2 on the VLT, A1FOSC on the Nordic Optical 
Telescope, GMOS-N and GMOS-S on the Gemini North 
and South, respectively, and finally LRIS on the Keck 
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Telescope. Depending on the grating or grism chosen 
for the observations, these instruments provided spectra 
with resolutions between 315 and 2140. We exclude after- 
glow spectra with uncertain redshift determination due 
to the lack of absorption lines in their spectra, and where 
the redshifts were determined from photometric redshift 
methods from their host galaxies. We also exclude af- 
terglow spectra where only upper limits on the redshifts 
have been established. In total 60 afterglow spectra pass 
these criteria and were used to create the composite spec- 
trum. Their redshifts lie in the range 0.35 < z < 6.7, 
with a mean of z = 2.22. 

Given the variety of instruments used, time after the 
trigger for the spectroscopic observations following the 
burst, and intrinsic magnitude of the afterglow, the S/N 
of the individual spectra vary between a few and ^200 
for the burst with the highest signal (GRB 060729 at z = 
0.5428). 

For the z > 2 GRBs, F09 determine their neu- 
tral hydrogen column densities to lie in the range 
17.0 < log(Af(Hi)/cm~ 2 ) < 22.7, and most of them 
are in the DLA regime. Metallicities have been deter- 
mined for a subsample of afterglow s which give typical 
values between 1% and 30% solar dWatson et alj|2006t 
Fvnbo et all 120061: iProchaska et al.l l2007ri iFox et al.l 
2008tlEhasd6ttir et al.ll2009l : lProchaska et al.ll2009h . The 
associated metal absorption lines are typically very 
strong and are certain to be associated with the host 
galaxy itself. 

Although the observational biases for the sample as a 
whole are well known (F09), the effects on the composite 
spectrum can be severe. For example, dust obscuration 
will render the afterglow fainter, and the S/N will typi- 
cally be lower. Bursts which are observed early, have low 
extinction, lie at lower redshifts, occur in a certain envi- 
ronment, or also possibly in an environment with a spe- 
cific chemical composition at a fortuitous distance from 
the burst, are more likely to have a bright afterglow. Be- 
cause of the weighting scheme for creating the composite 
spectrum these bursts will dominate over the remain- 
ing ones. Alternative combination schemes are found to 
produce a spectrum with a worse S/N and one with a 
smaller useful wavelength range, but general properties 
are found to be consistent, so this paper focuses mainly 
on the weighted combination scheme described below. 

3.2. Composite spectrum 

Each individual spectrum was normalized with a cubic 
spline function. The continuum level was estimated by 
fixing spline points in regions that were not contaminated 
with any absorption lines, either intrinsic to the GRB or 
from intervening systems as listed in the tables in F09, 
or from the atmosphere. The regions for the continuum 
was estimated by eye, and special care was taken in the 
wavelength regions around and blueward of the Lya line. 
This manual procedure was estimated to be most appro- 
priate since the individual spectra have widely different 
quality and spectral range, so that an automatic proce- 
dure would be very complicated. Since the normalization 
was done by eye, the errors introduced on a pixel-by-pixel 
scale are not known, so we cannot propagate the error 
from the normalization procedure throughout. With it- 
erations of the normalization, we estimate from the spec- 
tra with the best signal that we can achieve an accuracy 



better than 1% for the normalization. Whenever pos- 
sible, we used the flux calibrated afterglow spectra to 
derive the normalized spectra, because these should fol- 
low a pure power law. For 15 of the 60 spectra, we did 
not have a flux calibrated spectrum, while other flux cali- 
brated spectra had artificial features introduced by errors 
in the flux calibration from the spectrophotometric star 
as mentioned in F09. These features were removed too 
by the normalization. 

The associated error spectra were normalized as well 
using the same spline function. GRB afterglows can have 
different spectral slopes, so the normalization was a nec- 
essary step before the composite spectrum is computed. 
Any specific signature of the GRB environment such as 
dust extinction, and the 2175 A bump that appeared 
for example in the GRB 0708 02 afterglow (Kriihlc r et al.l 
I2008t lEhasdottir et al.ll2009H . is eliminated during this 
process. Here we do not aim to analyze the composite, 
line of sight averaged, extinction curve. 

Before combining the spectra, we masked absorption 
features that were unrelated to the GRB environments. 
To avoid contamination by absorption lines from inter- 
vening absorption systems, we masked each of the de- 
tected lines arising in intervening systems which were 
not associated with the GRB (see Tables 5-72 in F09). 
We masked fringes present at red wavelengths in Al- 
FOSC spectra, and fringes present at wavelengths larger 
than 7000 A in spectra obtained with the blue sensitive 
FORS1 detector installed in April 2007. We also masked 
atmospheric features such as the strong telluric bands 
around 6900, and 7600 A, but also weaker telluric bands 
around 6300, 7200, 8200, and 9000 A were masked; if not, 
they would introduce artifacts and spurious absorption 
lines redward of 3000 A in the rest frame. The masking 
was done by assigning a large value (10 6 ) to the error 
spectrum at the corresponding wavelengths. 

The typical dispersion obtained with the FORS V300 
grating is 3.25 A pixel -1 and the average redshift of all 
the bursts which contribute to the composite is z = 2.22. 
The average redshift therefore matches the dispersion of 
the low resolution perfectly A res t = A b s /(1 + z), so the 
optimal choice is to shift each spectrum to z = with a 
dispersion of 1 A pixel -1 . We investigated the outcome 
of the composite spectrum by choosing other resampling 
sizes for the dispersion, e.g., 0.5 and 0.8 A pixel -1 , how- 
ever, the absorption lines were not better resolved and 
the S /N decreased correspondingly. We adopted the red- 
shift for each GRB from neutral species associated with 
the GRB environment as listed in the tables in F09. 

The spectra were combined using the weighted average 

/A.comb = E^Ka)/ £ Vof,A 

i i 

where the flux / at each wavelength A was weighted by its 
error er, and summed over all spectra, i. The associated 
error spectrum was also calculated: 

4A,co m b = (E 1 /^,A)- 1 - (2) 

i 

The resulting composite spectrum and its error spec- 
trum is shown in Fig. [TJ and Table [1] lists the relative 
fluxes and uncertainties. The median value of the er- 
ror spectrum is 0.0065 from 1000-2000 A, corresponding 
to a median S/N = 150. At wavelengths 2000-3500A 
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the error is smaller, and the S/N = 250-300 per A, as 
illustrated in the middle sub-panels in Fig. [TJ At wave- 
lengths larger than 3500 A the S/N depends on the ex- 
act region, where only few, low signal data contribute, 
and where some regions have lower signal because the 
input spectra receive a low weight in regions of strong 
telluric absorption lines. Only one low redshift afterglow 
(GRB 061021 at z = 0.3463) contributes to the compos- 
ite beyond 5800 A, and here the noise is significantly 
larger than in other spectral bins. 

Apart from the regular error spectrum there are other 
effects that add to the noise. The composite spectrum 
is affected by uncertainties from the normalization of 
the individual spectra. If this is not done correctly, we 
can artificially introduce absorption lines in the compos- 
ite. At the very high S/N level, it is necessary that 
the normalization is done to a level better than 0.5- 
1%, otherwise features which resemble absorption lines 
will be introduced to the composite spectrum artificially 
and the uncertainty will introduce systematic errors for 
the equivalent widths. The normalization is particu- 
larly important for the individual spectra with good S/N 
which dominate the composites, however, when the S/N 
is good, the continuum flux level is also easier to identify 
visually. 

As afterglow spectra show a range of abundances and 
equivalent widths for individual absorption lines, it is 
interesting to determine the spreads of values with wave- 
length. We examined the variance of values from the 
sub-set of the 60 spectra that went into each wavelength 
bin. This variance spectrum essentially shows a flat be- 
havior with a value around 0.05-0.1 from 1000-4000A, 
while redward of 4000 A it has a larger scatter of val- 
ues between and 0.5. The reason why the variance 
spectrum does not have a larger value at the location 
of the metal absorption lines is likely because many of 
the individual spectra have low S/N, so they will domi- 
nate the spreads of values. One exception exists at the 
Lya line around 1216-1230A, where the variance is larger 
(0.3-0.4) because of the large range of hydrogen column 
densities probed by individual GRBs. 

A few individual high S/N spectra will necessarily 
dominate the composite spectrum. In order to check 
for systematic effects we also created a simple median 
spectrum, again masking out intervening absorption sys- 
tems as well as atmospheric absorption lines. The me- 
dian spectrum has a lower S/N ~ 50, estimated by the 
standard deviation in regions free of absorption lines, and 
a useful spectral range between 1000-2700 A. 

In addition to this, we also investigated the output 
where the spectra were multiplied with each other and 
divided by the number of systems. In practical terms, the 
spectra were added in logarithmic space, while omitting 
the masked pixels. 

l0g(/Axomb) = S i l0g(/ 4 ,A)/A r A, (3) 

where N\ is the number of spectra which contributes 
to each wavelength bin. Without the normalization, the 
composite spectrum will be completely dominated by ab- 
sorption lines, and have a very low output S/N. 

Physically, this method corresponds to placing the ab- 
sorbing systems one after the other along the line of sight, 
and then taking the average. In this case we avoid the 



systematic effects from weighting the individual spectra, 
however, we are again limited by the resulting S/N which 
is also <50. This implies that many of the faintest ab- 
sorption lines are lost in the noise. 

4. RESULTS 

Unless specifically mentioned, the following sections in- 
volves the weighted combination of the 60 afterglows, 
because we found that the general results did not de- 
pend on which of the three combination methods we 
used. The largest amount of information was therefore 
obtained with the highest S/N spectrum with the largest 
rest-frame wavelength coverage. Some results were also 
derived by creating sub samples from the full sample of 
bursts. 

4.1. Absorption line list 

Compared to individual spectra, we see significantly 
more absorption lines in the composite. Even the high 
resolution afterglow spectra obtained, e.g., with UVES or 
HIRES, are not able to detect some of the fainter lines 
at a high significance; although high reso lution spectra 
are b etter suited to investigate lines (e.g.. iLedoux et al\ 
120091 ) which are blended at the lower resolution. 

Table [2] lists all the lines and rest-frame equivalent 
widths (W r ) detected at above the 3c uncertainty level. 
Among all the lines identified, we checked if there are 
any missing lines that should be detectable. From the 
singly ionized species which are the dominant ionization 
state, we only find three non-detections, all of which are 
expected to be very weak. The non-detection of Nin 
A1502, Con A1466, and Cm A2040 are therefore not 
surprising. Some afterglow spectra have unidentified ab- 
sorption lines, as several ones in F09 do. After a care- 
ful masking of telluric features in the individual spec- 
tra, the composite spectrum does not appear to have 
any unidentified lines. Other very weak absorption lines 
which were uniquely detected in the GRB 080607 spec- 
trum, such as Oi A1355, Mel A1 683, Con A1574 and 
Gen A1602 ()Prochaska et al.ll2009|) . are not recovered in 
the composite spectrum. 

Identified lines which have rest-frame wavelengths sep- 
arated by 2 — 5 A were deblended by fitting Gaussian 
functions to the blends of lines using the splot routine 
in I RAF, keeping the wavelengths of the lines fixed, and 
fitting a single absorption line width. This routine also 
estimated the uncertainties using a representative value 
from the composite uncertainty spectrum around the line 
blend to represent the pixel statistics. Lines closer to- 
gether than ^2 A were not deblended and here the to- 
tal W r is listed. Reference rest-frame wavelen gths for 
the tr ansitions were found in the atomic table in lMortonl 
(2003). The W r uncertainties do not include the uncer- 
tainty for the placement of the continuum. 

Since some spectra with high S/N dominate the com- 
posite, we compared with the median composite and the 
multiplied spectra. Although their overall S/N are lower, 
we still find the same absorption lines in the three com- 
posites and the W r of the strong absorption lines are 
consistent to within the errors. 

To check further whether some of the weaker lines are 
randomly caused by errors in the normalization of the 
individual spectra, we created two independent compos- 
ites. The spectra were sorted according to their S/N, and 
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Fig. 1. — Each sub-spectrum consists of three parts. Top: The composite spectrum and its associated error spectrum (red line) which 
typically has a value around 0.005. Middle: The corresponding S/N per pixel, and the Bottom: shows the number of afterglow spectra 
that are used in the composite per spectral pixel. All detected absorption lines from Tables l2l and l3l are indicated. Note the change of scale 
on the y-axis, which implies that the error spectrum is not always visible. 
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two subsets of composite spectra with the same S/N were 
created. Since the GRBs have different redshifts, the best 
method would be to make two composites with a simi- 
lar S/N ratio throughout the spectral range, and with a 
similar redshift distribution so that there are an equal 
number of objects contributing to each spectral pixel. 
This is not feasible in practice, so the spectra were sim- 
ply combined according to their S/N ratio, and the two 
spectra have almost the same average redshift z = 2.4 
and z = 2.0, respectively. We checked whether the ab- 
sorption lines found in the full composite spectrum are 
present in both of the sub composites. This is not always 
the noted in the comments to Table [21 We can- 

not rule out that the lines are not real based on this test; 
the absorption lines may potentially be present in very 
few lines of sight, but randomly one that was observed 
at a high S/N. Additionally, the presence or absence of 
the absorption lines can be due to temporal variations as 
discussed in Sect. 14.41 In particular, some lines appear 
quite strong in one composite, but are absent in the other 
composite spectrum. 

Due to the wide range of iV(Hi) column densities, the 
Lya absorption line does not have a purely damped line 
profile, and residual emission can be seen in the Lya 
trough. Nevertheless, the average value of logTV(Hi) = 
22.0 provides a good Voigt profile fit to the red wing of 
the line in the composite spectrum. At this high column 
density, the red wing of the absorption line extends to 
1400 A which is consistent with the theoretical line pro- 
file. It also affects the continuum estimate of the metal 
absorption lines in this wavelength range. 

Blueward of the DLA, a mixture of metal absorp- 
tion lines and intervening hydrogen absorption lines are 
present. Hence, measurements of the equivalent widths 
of lines in this region are likely to be contaminated by 
intervening Hi clouds. Identified lines with wavelengths 
< 1215 A are listed in Table El 

As a consequence of the low resolution of the spectra, 
we see no evidence for any velocity shift of the absorp- 
tion lines relative to each other or multiple components of 
individual lines as reported in studies of individual after- 
glows observed wit h higher resolution spectroscopy (e.g. , 
Moller et al.1 12001 iFioreeTall [20051: iChen et all 120051 : 
Prochaska et al . 2008b; Fox et al. 200|). 



4.2. Lya forest transmission 

At z ~ 2, the intergalactic medium (IGM) expresses 
itself in quasar spectra as a series of H i Lya transi- 
tions commonly termed the Lya forest (|Rauchl 11998 - ). 
In a single spectrum, the Lya forest appears as a series 
of discrete absorption lines. If one considers an ensem- 
ble of sightlines, however, the stochastic nature of ab- 
sorption produces an average, effective opacity T e g. Us- 
ing hi gh S /N spectra of bright quasars iDairAglio et al.l 
(2008) fits the transmission in the Lya forest as a func- 
tion of redshift by the equation T = e~ r ° ff where r c ff = 
10-2- 21±0 - 09 (l + z )3.04±o.u In contrast to quasars w hi c h 

ionize their surroundings and create a proximity region 
with higher transmission nearby, GRBs and their hosts 
are expected to have relatively less impact on their sur- 
roundings, and the proximity effect will not be as pro- 
nounced. Therefore, GRBs are excellent probes of the 
transmission in the Lya forest, in particular, in the view 
that they can be observed at extremely high redshifts 



where the current record holder is GRB 090 423 at z = 8.2 
(jTanvir et al.ll2009l : ISalvaterra et al.ll2009f ). 

At the average redshift of z = 2.22 for the composite 
spectrum, the expected transmission in the Lya forest 
is T = 0.8 ±0.1. In the composite spectrum, many of 
the absorption lines in the Lya forest are associated with 
the ISM in the GRB host galaxy (see Table © . Without 
correcting for these lines, we find that the transmission 
is 0.6, but if we avoid the associated lines, the regions of 
the largest transmission where no metal absorption lines 
are expected suggests T ~ 0.85, as illustrated in the first 
two panels in Fig. [T] This is in excellent agreement with 
the prediction from the quasar transmission. Thus, we 
find no overall difference in the Lya forest transmission 
towards GRBs and quasars. To analyze the forest and 
its evolution in further detail, one needs to look at higher 
spectral resolution data. 

4.3. Lyman limit emission 

The existence of the Lya forest demonstrates that the 
IGM is predominantly ionized at z < 6. While quasars 
are known sources of ionizing radiation, they may not 
dominate t he extragalactic UV backgrou nd, especially at 
z > 3 (e.g.. iFaucher-Giguere et al"1l2009[ ). Instead, star- 
forming galaxies may contribute the majority of ionizing 
photons, but this requires that they have a non-negligible 
escape fraction / csc . Since GRB hosts are potentially 
unbiased tracers of star formation they have been used 
as probes of the ionizing radiation escaping blueward of 
the Lyman limit into the IGM. The estimated escape 
fraction from the Hi column density distribution func- 
tion for GRB afte rglows at z > 2 is (f csc ) = 0.02±0.02 
(|Chen et al.ll2007l ; F09). This value is consistent within 
the la error with numeri cal simulations, whi ch predict 
an escape fraction of 1% ([Pontzen et al.ll2010[) . 

Blueward of the Lyman limit in the composite spec- 
trum, we measure a relative flux level of 0.0060±0.0050 
in the wavelength range from 800 to 900 A as illustrated 
in the small inserted panel in Fig. [1] At these rest-frame 
wavelengths, only 11 afterglow spectra contribute to the 
composite, and some only cover redward of 880 A and 
have different S/N ratios in this wavelength range as 
listed in Table [11 Most afterglows at z > 2 have H I 
column densities in the DLA regime, where the escape 
of ionizing radiation is effectively blocked. Indeed, the 
average H I column density from the bursts in Table [H 
is 10 22 1±01 cm -2 , and at this high value the radiation 
blueward of the Lyman limit is completely absorbed. 

The average redshift for those afterglow which con- 
tribute blueward of the Lyman limit is z = 3.66. At 
this redshif t the average transmis sion in the Lya forest 
is T = 0.52 (IDairAglio et al.l l2008D. and in the Ly/3 forest 
it is T = 0.57 (|Madau|[T995). Including these two trans- 
mission corrections, the measured relative flux blueward 
of the Lyman limit is / = 0.020 ± 0.017, which is consis- 
tent with the values derived from the H I column density 
distribution. 

Since clouds with column densities in the DLA regime 
effectively absorb all the Lyman limit photons, a signif- 
icant escape fraction is only expected to be present for 
low column density systems. A few afterglows have re- 
vealed column densities well below the DLA limit, and 
GRB 050908 with log (iV(H i) /cm~ 2 )=17.6±0.1 has a sig- 
nificant level of ionizing radiation detected beyond the 
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Lyman limit (/ eS c = 13±4%), and appears to be a unique 
object among the whole sample. Because of the relatively 
low S /N for the GRB 050908 spectrum, eliminating it 
from the total composite changes flux blueward of the 
Lyman limit very little: / = 0.018 ± 0.017. A couple of 
other systems observed with low resolution spectroscopy 
(GRB 060124, GRB 070411, and GRB 071020) have rela- 
tively low Hi column densities, but their spectral ranges 
do not reach the Lyman limit. More studies which cover 
the very blue end of the spectral range are needed to de- 
termine the effects of the ionizing radiation from GRB 
host galaxies, and their impact on the extragalactic UV 
background radiation field. 

4.4. Temporal variations 

Our composite spectrum ignores the case of time 
variations in the fine structure levels. Several spec- 
tra show that the absorption lines strengths can 
change with time due to the temporal variation 
of the ionizing radiation from the central explosion 
(iDessauges-Zavadskv et ail l2006t IVreeswijk et al.l 120071 : 
iD'Elia et al.ll2009at iLedoux et al.ll2009j ). 

The spectra included in this investigation are obtained 
at one particular time (Table 2 in F09) at an average of 3 
hours after the burst in the rest frame for the composite 
spectrum. Since the spectra contribute to the composite 
weighted by their signal, the effective rest-frame observ- 
ing time can be calculated by weighting according to a 
representative value for the S/N of each spectrum. This 
gives an effective observing time of i ros t = 2.78 hours 
after the burst in the rest frame for the composite spec- 
trum. We investigate the cumulative effect of the observ- 
ing time on the spectra, and look for temporal variations 
in absorption lines. We create two subsamples separated 
according to the time of the start of the spectroscopic ob- 
servations after the burst trigger. The division takes into 
account the redshift of the GRB by shifting the start of 
the observing time to the rest frame. These sub-samples 
of composite spectra include 30 afterglow spectra ob- 
tained within 1.9 hours from the trigger (in the rest frame 
of the burst) which has an average redshift of z = 2.4, 
and 30 spectra obtained between 1.9 and 14 hours, which 
have an average redshift z = 2.0. Both composite spec- 
tra have a S/N between 100 and 200 in the wavelength 
range 1500-3500 A. The magnitudes of the afterglows at 
the time of the trigger vary by ~ 4 mag in both the early 
and late-time samples, while the early sample is on av- 
erage 1.4 mag brighter than the late sample (See Fig. 2 
in F09). 

A clear difference in the strength of the metal absorp- 
tion lines is seen. Strong variations in the iron fine- 
structure lines are observed from the spectrum obtained 
at an average of 0.4 hours after the burst to that obtained 
5 hours after the burst. However, also the species Fell, 
Aim, Mgn, Znii, Cm, Nin which trace the neutral gas 
seem to decrease in strength. In particular, the Fe n W r 
decrease by a factor of 2-4 as illustrated in Fig. [2] Be- 
tween the two composites, all the absorption lines ap- 
pear to decrease in strength by a factor of ~2. Other 
lines like C II, and Si II and the high ionization lines Si iv 
and C iv remain constant within 10-20% as illustrated 
in Fig. [3] The disappearing CO lines in Fig. [3] could also 
be ascribed to a temporal variation, but as explained in 
Sect. l5.3l the CO lines are only detected in the afterglows 



from GRB 080607 and possibly also GRB 060210. Hence, 
it is not certain whether this is a real effect or just a co- 
incidence from the strong CO lines in one GRB sight 
line. It is curious that GRB 080607 also has the earliest 
rest-frame observation time in the entire sample. 

The simple median spectrum and the multiplied spec- 
trum composites do not suffer from the possible biasing 
effects of weighting with noise level. Nevertheless, they 
also reveal changes in the absorption line strengths (both 
ground and excited states) between the early and the late 
time samples. 

To date, no time series of individual afterglows have re- 
vealed any significant variations in the non-exited states. 
In particular, a time series of high resolution spectra 
of GRB 071031 showed constant high ionization lines 
(|Fox et al.l I2008T ) . Therefore, we investigate whether 
other effects are at play in the composite spectra, and 
find the temporal changes of absorption line strengths to 
be a random effect caused by some afterglows with very 
strong absorption lines. In order to conclude this, we 
investigate Wjr sn and WfJ ou from individual systems as 
a function of the starting time for the spectroscopic ob- 
servations as well as a function of redshifts as shown in 
Fig. SJ First, there are no systematic changes of the W r 
with redshift. Second, there is a tendency that earlier 
observations have larger W r , but including the measured 
W r uncertainties any statistically significant correlation 
disappear. Both the early- and the late-time observations 
have spectra with high S/N ratio. Taking the W r at face 
value in addition to the fact that the earlier observations 
have brighter afterglows and therefore give more weight 
to the composite, this mimics the behavior of a tempo- 
ral variation of the composite spectrum. We therefore 
suspect that this effect is not physically real. A similar 
analysis of the C iv and Si iv W r measurements and their 
S/N show a large scatter of values at all times, so the ap- 
parent constant value in the composites is also likely to 
be a random effect. 

Apart from the composition of the external medium, 
the fine-structure line variations also depend on the in- 
put UV flux from the GRB. We investigated the effect 
by analyzing subsets of composites by separating accord- 
ing to the spectral slope /3, derived from the power-law 
dependence of the flux f v oc . Values of (3 are listed 
for half of the bursts in iKann et alj (|2010T ). while the 
remaining ones are obtained by fitting a power law to 
the flux calibrated spectra without including additional 
extinction. Including a set of possible extinction laws 
allows the spectral index a large range of values within 
the uncertainties. Here, we use the face values without 
including the uncertainties and create two subsamples 
separated at /3 = 1. Comparing the composite spectrum 
with lower fj values (brighter UV than optical flux in fx 
units) with that obtained from afterglows with j3 > 1 , we 
find that both have strong fine-structure lines, with no 
clear trend that one has systematically much stronger 
lines than the other. A more careful treatment of the 
spectral slope dependence would need to include system- 
atic analyses including the uncertainties and the intrinsic 
extinction to derive the intrinsic UV flux from the after- 
glows. Such an analysis is beyond the scope of this paper. 

4.5. Effect of GRB energy release 
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Fig. 3. — The composites show smaller variations within 20% for the high ionization lines C IV AA1548,f 550 and S IV AA1393,1402. Note 
the disappearing CO lines and weakening of the Si II* line. 



A trend between the GRB isotropic en ergy release 
and t he afterglow brightness has been found ()Kann et al.l 
2010), albeit with a very large scatter possibly due to 
variations in the circumburst density. Because of this 
possible connection, we explore the effects for the com- 
posite spectrum by the GRB energy release. GRBs and 
their afterglows are thought to be collimated with a typi- 
cal beaming angle between 1-25° (jFrail et al.l l2001). Out 
of the 60 afterglows investigated here, 42 bursts have a 
measured bolometric isotropic energy (|Kann eTal1l2O"l0L 
and references therein). These 42 bursts have energies 
between 51.38 < logi?i SOi boi(erg) < 54.49. 

We investigated the effects for the composite spec- 
trum by dividing the sample into two samples accord- 
ing to the estimated bolometric, isotropic energy re- 
lease from the prompt emission. Each subsample con- 



tains 21 GRBs. The energy at which the sample is di- 
vided is determined by statistical sampling rather than 
a physically real property. The low and high energy 
composites have 51.38 < log i?i SO ,boi (erg) < 52.88 and 
52.88 < log -Ei so ,boi (erg) < 54.49, respectively. Overall, 
there is no significant difference between the absorption 
line widths for either the high or low-ionization lines, and 
fine-structure lines are present in both the high-energy 
composite and low-energy composite. However, redward 
of 2200 A, all lines (Fe, Mn, Mg) in the low-energy com- 
posite are a factor of ~4 weaker. The reason for this is 
most likely due to the fact that a couple of high-energy 
GRB afterglows at z ~ 1 (GRB 061121 and GRB 080413) 
have rather strong Fe absorption lines and dominate the 
absorption lines at rest-frame wavelengths larger than 
2200 A. Neglecting these two GRBs, the high- and low- 
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energy composite spectra are similar, and the variations 
between the lines widths are 11%±22% with no system- 
atic trends of which of the low- or high-energy compos- 
ites produce the strongest lines. We conclude that the 
prompt GRB energy release itself has no effect on the 
absorption lines. 

4.6. Dark bursts 

For some bursts, afterglows are not detected, but the 
success of detecting the afterglows depend strongly on 
the afterglow search strategies, i.e., rapid searches which 
are able to probe faint magnitudes. Some bursts may be 
dark, because they lie at high redshifts and are therefore 
not detected in optical follow-ups. Others may be ob- 
scured by dust. A better method to define intrinsically 
dark bursts relies on the X-ray flux versus the simulta- 
neous optical aftergl ow magnitude. From th is optical to 
X-ray spectral index, Uakobsson et al.l (|2004h define dark 
busts as those having /3ox < 0.5. 

In the sample of 60 afterglows, 53 have measured /?ox 
(Table 73 in F09), and 10 of these satisfy the criterion 
for being a dark burst. We created composite spectra for 
those with /3ox < 0.5 and /?ox > 0.5, but because of the 
large difference in the number of spectra that go into the 
composite, the S/N is also unequal. The composite spec- 



trum from the bright burst sample with /?ox > 0.5 has an 
S/N between 100 and 200 from 1000 to 3000 A, while the 
composite from the dark burst sample with /3ox < 0.5 
has a S /N = 25-40 in the same wavelength region. Still, 
the dark burst composite has stronger absorption lines 
by a factor of 2-3 over this spectral range, while Hi col- 
umn density decreases by 0.2 dex from the bright to the 
dark composite. A representative wavelength region is 
shown in Fig. [5] 

While the optical to spectral index separation at /?ox = 
0.5 is based on a physical reasoning, we may still won- 
der if artifacts are included from comparing two spectra 
with very different S/N ratios. Choosing instead a dif- 
ferent separation of /3ox for the subsamples may involve 
bursts with slopes that are dependent on the particu- 
lar positions of the spectral breaks between the optical 
and X-ray bands. Although the separation is not aris- 
ing from a physical explanation, the low-/3ox sample will 
still include the 'darkest' bursts. Dividing the full sam- 
ple into two equal size parts instead with a separation 
in the spectral index of /?ox = 0.75 such that the S/N 
is comparable, consistently shows that the darker bursts 
have stronger absorption lines. The detection of stronger 
metal absorption lines towards dark bursts is fully con- 
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sistcnt with the idea that these reside in more metal-rich, 
dustier ISM relative to the bright counterparts. This is 
already established from analysis of the pro perties of the 
host galaxies from broad band photometry (|Perlev et al.l 
2009) and from X-ray absorption (F09). 



4.7. Afterglow optical luminosities 

The optical luminosities of the afterglows depend on a 
large number of parameters, including the density of the 
external medium. Although the medium which causes 
the absorption lines in afterglow spectra is generally de- 
termined to lie at distances greater than 1 pc, we may 
still investigate if the composite spectrum shows any 
changes related to the optical luminosities. Potentially, 
some absorption lines may arise closer to the GRB ex- 
plosion than otherwise determined, or there may be an 
anti-correlation between a more dusty, metal rich envi- 
ronments and the optical luminosity of the afterglows. 

Having investigated the effect of the spectral slopes, we 
include an analysis of any systematic effects of the after- 
glow optical br ightness. Analyzin g the light curves in a 
systematic wav. lKann et al.l (|2010f) determine the observ- 
able, extinction corrected magnitudes of a large number 
of afterglows at t = 1 day from the explosion in the rest 
frame of the bursts. The Rc magnitudes are listed for 
39 of the 60 bursts in the composite sample. We divide 
the spectra into two almost equal sized subsamples sepa- 
rated at Rc = 19.55 mag. The two resulting composites 
have the same S/N ratio. 

The optically faint subsample has stronger Fen lines 
redward of 2200 A compared to the bright subsample. 
While the high ionization lines in the two samples are of 
the same strength, the low-ionization ones are about 50% 
stronger in the faint sample. As for the isotropic energy 
release case, two bursts (GRB 061121 and GRB 080413) 
are among the optically faint ones 1 day after the explo- 
sion. At the time when the spectra were obtained they 
were among the brighter ones. Eliminating these two 
GRB afterglows from the faint sample, the composite 
spectra have similar absorption line strengths. 



5. METAL ABSORPTION LINES 

In this section, we compare the individual features 
in the composite spectrum with other categories of ab- 
sorption line systems with the purpose of investigating 
whether GRB environments are special. This is interest- 
ing in the context of constraining the nature of the GRB 
progenitor through any peculiarities in the surrounding 
medium. 

5.1. Strong Can and Mgn absorbers 

The composite spectrum displays a large number of ab- 
sorption lines. To compare the composite with other high 
redshift galaxy spectra, we arc limited to the stronger 
transitions found in large samples in low spectral resolu- 
tion data, because these are the ones most easily detected 
in faint galaxy spectra. 

The Can absorption lines are very strong in the com- 
posite with W r Ca " A3934 = 0.77 A, which places the av- 
erage GRB site in the category of the dustiest absorbers 
(jZvch et alJl2009h compared to other intervening absorp- 
tions systems selected by their neutral hydrogen column 
densities. In intervening quasar absorbers at lower red- 
shifts, strong Can lines b elong to DLAs w ith a color ex- 
cess of Eb-v > 0.1 mag (jWild et al.ll2006| ). The reverse 
is not true; only 1 0% of DLAs have strong Can lines 
(|Nestor et al.l [2()()M. and DLA selected absorbers have 
a n order of magnitude les s reddening {Eb-v < 0.02 mag 
in lMurphv fc Liskdl200l . 

Also the Mgn absorption lines fall among strong, 
but not extreme a bsorbers detected in quasar spectra 
(jNestor et all 12005( 1 . Its W* 2797 = 1.7 A is comm on for 
systems associated with DLAs (Ra o et ail [2006) . The 
Mgn W r is smaller than the W r = 3.8 A found in a com- 
posite spectrum of local starburst galaxies (C. Tremonti, 
private communication) as well as ultra strong (W r > 3.0 
A) Mgll absorbers at z = 0.7 dNestor et al.ll2010l) , whic h 
possibly trace starburst galaxies (|Rubin et al.1 l2010f l. 
The widths of Mgn lines in quasar spectra which arise 
in intervening galaxies at 0.3 < z < 1 also depend on th e 
galaxy's rotation and inclination ([Kacprzak et al.l f2010). 
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Since GRB hosts have low stellar masses ( Savaglio ct al. 
2009), star burst driven galactic winds would easily es- 
cape their potential well. Strong lines with large W r such 
as Sin A1526 have very large optical depths such that 
the line widths measure the velocity spread of individ- 
ual compon ents in DLA systems b oth towards quasars 
and GRBs (|Prochaska et al.ll2008al) . Large widths can 
therefore be seen as proxies of galaxy wind strengths, 
which would indicate that GRB hosts have weaker winds 
relative to galaxies with higher absolute star formation 
rates (SFRs). Higher resolution spectra are needed to 
explore the effects of galaxy winds in GRB hosts in or- 
der to quantify the effect the host mass and SFRs on the 
absorption line widths of the very strong lines. 

5.2. Dusty afterglow spectra 

There is currently little evidence either for a strong ex- 
tinction from the af terglows themselves with typical val- 
ues o f E B -v < 0.3 (jGalama fc Wijersll200H iKann et al.1 
2006), while a few individual bursts suffer from more 
extin ction like GRB 07030 6 which had A v = 5.5 ± 0.6 
mag (|Jaunsen et all 120081 ) and GRB 070802 which had 
Ay = 0.8 — 1.5 mag and sho wed a pronounced 2175 A 
bump (|Eh'asd6ttir et al.ll2009j) . These two GRB spectra 
have relatively low S/N, so they do not affect the com- 
posite significantly. Recently, a very high extinction of 
A v > 12 was suggest ed for the dark burst GRB 0904 17B 
(|Hoiland et al.l l2010h . However, such dusty environ- 
ments are rare in the sample studied here, where the 
extinction measur ed from the afterglows is generally low 
(|Kann et al.ll2010ft . 

In the integrated spectral energy distribution (SED) 
of the host galaxies there is evi dence for relatively little 
reddening around En_v ~ 0.1 dChristensen et al.l [20041 : 
ISavaglio et aT1l2009l; iLevesaue et al.ll2010D . although fits 
to SEDs which include the far-IR e mission find larger 
amou nts of extinction in GRB hosts (|MichaIowski et all 

2008) . Analyses of individual afterglow spectral slopes 
generally only allow for small amoun ts of extinction 
in the burst region ([Kann et al.l I2010D . If no dust is 
present in the host galaxy, this would produce a brighter 
optical afterglow a nd introduce an observational bias 
(|Fvnbo et all [200ll V_Jbr_ some GRBs n o optical after- 
glow is detected ([jakobsso n et al.l [2004D . which is most 
likely due to dust obscuration in their host galaxies 
(|Perlev et al.l |2009L ; F09). However, for the compos- 
ite spectrum of bursts with detected optical afterglows, 
dust obscuration is not a strong effect. 

5.3. Molecular lines and diffuse interstellar bands 

We search the composite for absorption lines aris- 
ing from molecular lines, but only CO lines are found. 
These are very rare in GRB sight lines, and the com- 
posite is dominated by the strong CO lines observed in 
GRB 080607 and which are caused by excitation from the 
GRB (jProchaska et al.ll2009Hs"heffer et al.ll2009l) . A pos- 
sible, but weaker detection is that towards GRB 060210 
(F09). Eliminating the GRB 080607 afterglow while cre- 
ating the composite, the two detected CO lines in the 
composite vanish (see Fig. [3]). Regarding other molec- 
ular lines, H2 abso rption lines are not commonly found 
in GR B afterglows ()Tumlinson et al.ll2007t iLedoux et al.l 

2009) . The strong H2 lines observed in the line of sight 



to GRB 080607 are blended with too many other lines in 
the composite Lya forest, so we cannot include a more 
detailed study here. 

We also search for diffuse interstellar bands (DIBs) red- 
ward of 4000 A using a list of DIBs found in the Galaxy 
(jJenniskens fe Desert] Il994) . In particular, we look for 
the strongest bands at A4428, A5705, A5780, A5797, 
A6010 A6203, and A6283, which apart from seen in Galac- 
tic sight lines, are also found in nearby starburst galaxy 
spectra. Some of these lines have recently been observed 
in a high-resolution spectrum of the Type lb SN 2008D 
which exploded in NGC 2770 (jThone et al.ll2009D . show- 
ing that DIBs can be present in galaxies which form mas- 
sive stars. These lines are found to have typical strength 
of > 0.1 A, and their equivalent widths correlate with 
the color excess Eb-v and the hydrogen colum n den- 
sity (jHeckman fc Lehnertll2000l IWeltv et al.ll2006D . Even 
considering the very high iV(Hi) in the composite, we 
failed to get any significant detection although the spec- 
tral range and S/N ratio for the composite is sufficient 
to find some of the bands with W r ~ 0.1 A. IWeltv et al.l 
(|2006l) derive a best-fit relation between the 5780 A DIB 
strength and the reddening in Galactic sight lines. Using 
their fit logE B -v « -2.7+ 1.1 log Wi? (5780) /mA, the 
limit for W r corresponds to an upper limit of Eb-v < 0.3 
mag. A similar reddening limit would be inferred by a 
UV- shielded line of sight in the Small Magellanic Cloud 
(|Cox et al.M2007h . 

Relative to stars in the Galaxy, weaker DIBs are found 
in the less metal-ri ch medium of the Magellanic Clouds 
(|Weltv et all I2006D . A higher ionizing radiation field 
than that present in the Galaxy may be responsible for 
destroying t he carriers of the molecules that give rise 
to the DIBs (|Weltv et al.l 120061 : ICox et al.ll2007l ). Apart 
from detections of DIBs in the Ga laxy, DIBs have been 
found in some gas-rich galaxies. lLawton et al.l ()2008f ) 
find that strong DIBs are absent in six out of seven 
DLAs at z < 0.5, with the one detection belonging to 
the most metal rich syst em with the highest reddening 
(jJunkkarinen et aTll2004l ). Like the DLA systems, the av- 
erage GRB host ISM is also metal-poor and contain rel- 
atively little dust. Additionall y, GRB regions are k nown 
to have strong radiation fields ( Levesque et al.l2010T l . and 
therefore the DIB carriers are likely to be destroyed in 
the region near the GRBs. These arguments combined 
may be the reason for the non-detection of DIBs in the 
composite spectrum. 

The most commonly found DIBs fall in the rest frame 
visible region of the composite spectrum which has a sig- 
nificantly smaller S/N ratio, and where only few after- 
glows contribute. It would be more useful to look for 
DIBs in the UV part of the spectrum, which allows the 
detection of weaker lines. Some unidentified UV features 
with equivalent widths of a few to a few tens of mA 
possibly belonging to diffuse ban ds have been identified 
in spectra of Milky W ay stars (|Pwa fc Potta sch 1986: 
iDestree fc Snowj|2009fl . Of the few UV DIB candidates 
suggested, most lines lie too close to other metal absorp- 
tion lines in the low-resolution composite spectrum to 
allow their identifications. Two exceptions are UV fea- 
tures at 1384.3 A and 1490.1 A, which are not detected in 
the composite afterglow spectrum to a 3 a limit of 60 mA. 
The features observed in stellar spectra are even weaker 
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than this detection limit, so il UV DIBs are present in 
GRB spectra an even better S/N ratio is necessary to be 
able to detect them. 

5.4. High-ionization Lines 

In almost all low redshift starburst galaxies, the 
SiivA1400 and CivA1550 lines are quite strong with 
rest- frame equivalent widths >3A as seen in HST spec- 
tra of local starburst galaxies (C. Tremonti, private com- 
munication). In the afterglow composite, the absorption 
lines are weaker with W r —1 and 2 A, respectively. The 
transitions arise as a combination of stellar wind fea- 
tures and intrinsic ISM absorption. In starburst galaxies, 
the wind component is typically recognized from their P- 
Cygni profiles. Contrary to direct spectra of starbursts, 
the GRB afterglow spectra are not the direct emission 
associated with massive stars, so with the absence of a 
recognizable P-Cygni profile in the GRB afterglow spec- 
tra, a lower W r is expected. 

The C iv and Si iv lines typically extend over several 
hundreds of km in GRB spectra, and h ave several com- 
ponents some of whic h are saturated ([Prochaska et al.l 
l2008aHFox et al.ll2008l ; F09). In the low-resolution com- 
posite spectrum these lines have residual relative fluxes 
< 0.7, which suggests that they are composed of individ- 
ual saturated lines as explained in Sect. [7] 

The highly ionized gas could potentially trace the 
medium surrounding the GRB progenitors. Species 
like N V and O vi require a stronger ionizing radiation 
field than that produced by O stars, but which ca n 
be produced by a GRB event (|Prochaska et al.l l2008bj) . 
although the creation by collisional ioniza tion in the 
ISM cannot be ruled out (|Fox et al.l [2008h . The Nv 
AA1238,1242 doublet is distinguishable in some of the 
individual afterglow spectra at low resolution, while in 
high re solution spectra reveal that these lines a re ubiq- 
uitous (jProchaifeiran [20081 ^EeFaI][200l). These 
lines are not easily recognizable in the composite spec- 
trum, because they are blended with the red wing of the 
Lya line, in particular in spectra with high signals and 
very high Hi column densities. In the composite spec- 
trum the N V lines appear as a slight deviation from the 
broad wing profile of Lya. A more detailed analysis of 
these line requires a higher spectral resolution. 

6. COMPARISON WITH STARBURST GALAXY SPECTRA 

With a high S/N composite it is relevant to look for 
differences in the absorption line properties compared to 
what is known about other galaxies. With a spectral 
range covering from the UV to the rest-frame optical 
wavelengths, a comparison with local starburst galaxies 
is justified, but given the redshift difference it is more 
appropriate to compare with starburst galaxies at higher 
rcdshifts. Such a comparison is hampered by the faint- 
ness of the distant galaxies, but in a few cases a good 
signal from high redshift galaxies have been obtained. 

A low resolution composite spectrum from 811 Lyman 
break galaxies (LBGs), which has a measured S/N of 
~ 40 and a rest-frame wa velength coverage fro m 920 
to 2000 A is presented in IShaplev et all (|2003l) . The 
composite LBG spectrum has been created by averag- 
ing the spectra of galaxies at z > 3 with magnitudes 
R < 25.5, and the composite spectrum has a dispersion 
of 1 A pixel" 1 . With an average redshift z ~ 3, this is 



comparable to the GRB composite spectrum. Another 
example is the magnified signal from high redshift galax- 
ies lensed by foreground galaxies or clusters. The proto- 
type LBG spectrum obtained with Ke ck/ESI of MS 1512- 
cB58 (cB58 for short) at z = 2.7276 (jPettini et alJl2000l 
2002), which covers 1050—2800 A in the rest-frame and 
has an S/N~2-15, presents another valid comparison for 
the GRB composite. Finally, we also compare the GRB 
composite with a UV composite spectrum of UV bright 
galaxies at 1.3 < z < 2 from the G emini Deep Deep 
survey (GDDS) (jSavaglio et al.ll2004l ), which covers the 
wavelength range 2000-3000 A with an S/N of about 15. 

Fine-structure lines can arise in a medium which is be- 
ing pu mped by UV radiatio n, or in very high density 
regions (|Silva fe Vi cgas 2002]). Fell* fine-structure lines 
have been detected in several GRB sight lines but not 
in quasar absorption systems (|Chen et al.ll2005lh and lo- 
cally only in the dense wind from n Carinae ( Gull et ah! 
|2005| ). The UV pumping scenario and the time varia- 
tions of the absorption lines in time series observations 
of individual bursts indi cate a decrease of the ionizing 
radiation from the GRB (|Prochaska et al.ll2006l ). In this 
scenario, one would only expect the presence of Fell* 
fine-structure lines to be related to the environment in 
the GRB host galaxy, and generally not to be present in 
the ISM in normal starburst galaxies. 

Since l ong duration GRBs exp lode in star forming 
galaxies (|Christensen et all I2004D . we here investigate 
if the fine-structure lines could be present generally in 
starburst galaxies which have particularly large SFRs, 
or strong emission in the UV, or emission arising from 
dense environments. In a composite spectrum created 
from local starburst galaxies (C. Tremonti private com- 
munication) most of the metal absorption lines in the 
range 2150-3250 A identified in the GRB composite are 
also present (See Figure [6]). The strongest transition is 
that from Fen* A2396 which is clearly detected in the 
GRB composite, but due to the lower quality of the star- 
burst composite (S/N~30), we cannot claim a clear pres- 
ence of Fe II* fine structure lines. A better signal from a 
starburst galaxy is needed to establish if these lines are 
only present in the most extreme conditions in the region 
around a GRB explosion either by being pumped by the 
hard UV radiation, or by having a specially high density. 

In terms of tracing starburst galaxies, LBGs are 
brighter and more massive than GRB hosts at the same 
redshifts due to the selection of bright sources for the 
LBG samples, whereas GRBs trace star forming galax- 
ies irrespectively of their brightness. The same selection 
effects play a role for the galaxies in the GDDS, while 
cB58 has an intrinsic brightness of a typical L* galaxy. 
None of the comparison spectra show Fe n* lines; in the 
case of cB58 the lines around 1600, 2300, and 2600 A 
could have been detected if they had been there. The 
non-detection of fine-structure lines in starburst galaxies 
is perhaps not surprising given that Fe n* lines in after- 
glow spectra decay rapidly after the burst, and on that 
basis the lines are not expected to be present on the much 
longer timescales present for regions where massive stars 
are formed. 

With a few excep tions, the W r measured for cB58 in 
iPettini et al.l (|2002T ) have typical values about twice the 
strengths of those in the GRB composite. The exceptions 
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Fig. 6. — Section of the GRB afterglow composite spectrum (in black) compared to a starburst galaxy composite spectrum in red in 
the upper panel (from C. Tremonti). The dotted lines correspond to ltr deviations from the composite. While strong Fen* fine-structure 
transitions are present in the GRB composite, only the those from the ground state are present in the starburst composite. The lower 
panel shows the same region for the GDDS composite spectrum. 



are the Fen AA1608,1611 blend, Cm, Znii A2026 and 
Cm, Znn A2062 blends, and the Mnll A2577 line, all of 
which have roughly the same line strengths as in the GRB 
composite. In contrast, the LBG and GDDS spectra have 
in general line strengths in be tween the GRB com posite 
and cB58 (see also Table 2 in lSavaglio et al.ll2004D . 

LBGs and GDDS galaxies are relatively more massive, 
have higher SFRs, and have stronger outflows driven by 
massive stars and supernovae relative to the GRB hosts. 
The SFRs derived from the UV continuum flux relative 
to the GRB hosts lumi nosities are comparatively high 
(|Christensen et aT1l2004l ). In contrast the specific SFR 
(SFR per unit stellar mass) is found to be low compared 
to ot her high redshift galaxi e s (ICastro Ceron et al.ll2006l 
120101) , while ISavaglio et all ((20090 nn d the contrary. It 
is therefore still an open issue as to which degree the line 
strengths in GRB hosts and other high redshift galaxies 
are related to their specific SFRs. 

The C iv AA1548,1550 lines are stronger in all compari- 
son spectra. This line arises as a combination of ISM and 
stellar, outflow driven, absorption components. Since the 
GRB hosts are generally very faint, the emission from the 
host does not dominate, such that no P-Cygni profiles are 
expected in the GRB afterglow spectra. Other excep- 
tions are the Si n* lines which are found in absorption in 
the afterglow composite, while these lines are found in 
emission in the cB58 spectrum and the LBG composite. 

7. COLUMN DENSITIES 

The hydrogen column density for cB58 is about one 
order of magnitude below that in the GRB composite: 
log-ZV(Hi) = 20.85 resulting in a higher metallicity of 
40% solar (jPettini et alJl2002f >. Since the widths of the 
metal absorption lines in the GRB spectrum are gener- 
ally lower than in other high redshift galaxies, one is led 



to the question of what is the average metallicity in GRB 
lines of sight. However, this is a complex issue to ad- 
dress using the GRB composite spectrum, because each 
absorption line does not measure the same samples of 
individual GRB lines of sight. Any composite spectrum 
will carry along the scatter of properties of the individual 
objects (galaxies or GRB afterglows) that went into it in 
the first place, and a physical interpretation is hampered 
by this. Still, the general features in a composite spec- 
trum can be derived, as done from the LBG composite 
with its interstellar absorption lines (|Shaplev et al.|[2003l) 
and for the GDDS composite, where the trends of metal 
column densities, depletion pattern s and extinction wer e 
derived for galaxies at 1.4 < z < 2 (jSavag lio et al.l l2004). 
Compared to these studies, the GRB afterglow composite 
may present further complications, because the galaxies 
probed by the random sight line of a GRB through their 
ISMs are not as uniformly selected. 

Bearing this in mind we analyse the metal lines to de- 
rive column densities for a wide range of elements. Since 
the spectra are obtained at low spectral resolution, we 
use a standard single component curve-of-growth (COG) 
method to calculate the column densities. For GRB spec- 
tra which have absorption lines with multiple, saturated 
components spread over a velocity range of several 100 
km s _1 , the simple COG method underestimates the true 
metal column densities and the metallicities ca n be larger 
by u p to 1 dex in the worst case scenario (jProchaskal 
l200ll . 

In the case of Hi, its column density is determined 
instead from Voigt profile fitting to the red wing of the 
absorption line. In the COG fit, the best fit Doppler pa- 
rameter (b) is determined using ten individual Fell lines. 
Including all the uncertainties of the lines, a best fit yields 
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b = 44 ± 2 km s" 1 and l ogiV(Feli)=15.62 ±0.05. How- 
ever, as demonstrated in iProcha ska (200(|), the weaker 
lines in the linear part of the COG have underpredicted 
strengths, and the true column density may be larger 
than given by the best fit. Instead, we choose to tie the 
fit to the four weaker lines (Fen A1611, A1901, A2249, 
and A2260) as illustrated in Fig. [7] Since the widths 
of the stronger iron lines have been calculated by de- 
blending them from the Fen* fine-structure lines, some 
contamination may still occur, in particular for the Fe n 
A2382 and Fe II A2600 lines. The best fit is determined 
from x 2 minimization, which gives a Doppler parameter 
b = 34 ± 2 km s- 1 and log TV (Fen) = 15.70±0.07. All 
Fen lines follow a single curve well both in the linear 
part, while the points in the fiat part show some scatter 
around the best fit COG. To assess the amount of con- 
tamination by the fine-structure lines, we examine the 
composite spectrum created from the late-time sample 
in Sect. 14.41 We measure the W r for the 10 Fen lines 
and find the best fit with a smaller Doppler parameter: 
b = 20 ± 2 km s^ 1 and logiV(Feii) = 15.60±0.10. 

Since low resolution data indicate a higher b pa- 
rameter than when observed at high resol ution, be- 
cause it measures the spread of the lines (Prochaska 
2006) rather than the intrinsic velocity of the gas, it 
is relevant to compare with individual measurements 
of b from detailed high resolution spectra. Typical 
Doppler parameters measured in various GRB after- 
glow spectra are betwe e n 7 and 15 km s" 1 dFiore et al. 
20051 IChen etail [20071 [Thone et al.l l2007t iD'Elia et al. 
2009aibf ). while IVreeswiik et al.l (|2007t ) found b = 25 
km s _1 for GRB 060418. Hence, the reason for our find- 
ing of a higher b is partly due to blending with iron 
finestructure lines, while a minor remaining part could 
be caused by intrinsic line spreads of different compo- 
nents in a line. 

To derive column densities for other elements, we use 
the best-fit b parameter from the fit to the iron lines. 
Column densities for all the elements are listed in TableO 
We avoid using lines which are obviously blended with 
other elements as for example the O I A1302 and Si n 
A1304 blend. To derive the column density of Znii, we 
subtract the contribution of the CriiA2062 to the total 
width of the line. From fit of the other two Cr II lines to 
the COG we find that ^ r CrnA2062 = 0.17 A, whereas the 
Cr II A2026 line makes an insignificant contribution to the 
ZnllA2026 line. If the column density of Mgi is low as 
the best fit suggests in Sect [7l then the contribution of 
Mgi A2026 to the line blend is 0.05 A, which we subtract 
for the fit to the Zn n column density. 

We checked whether the composite spectrum is reli- 
able for the analysis of the metal column densities by 
comparing the fits with the spectrum derived from Eq. 
(3). Using the ten Fen lines, we get a consistent re- 
sult with the best fit of log iV (Fen) = 15.90±0.10 and 
b = 30 ± 2 km s _1 . With a lower S/N, this composite 
does not allow us to derive column densities of all the 
weak lines and rare elements. 

The elements are shown in Fig. [7] as different sym- 
bol shapes. This plot illustrates that in general, the 
absorption lines fit well to the Fell COG apart from 
at low values of log (TV/A (cm _1 ))<7, where the scat- 
ter may have a physical explanation. To mention a few 
elements in particular, the Tin lines are relatively weak, 



and the scatter for the six detected lines for this element 
in Fig. [7] is substantial. Also some of the Mg I lines are 
out-lying relative to the COG. The main reason is pos- 
sibly the sampling of the absorption lines by different 
GRBs. Since the Mgi lines lie in the range 1747-2852 A 
and Tin lines in the range 1906-3384 A, each absorp- 
tion line must sample a different set of GRBs. In these 
cases, we exclude points with log (NfX) <7, but warn 
that in reality the column densities for these elements 
could be ~1 dex higher than listed in Tabled The same 
tendency is observed for Fei lines. At the other end of 
the scale, Sin A1260 and Sin A1808 are blended with 
their respective fine-structure lines, and are not possible 
to deblend. As a consequence, the Si II lines in Fig. [7] 
scatter and we base the fit for the column density mostly 
on the isolated Sin A1526 line. Since the only singly ion- 
ized carbon line (C n A1534) detected in the composite is 
blended with Cn* A1335, we cannot calculate its column 
density accurately. Instead, we assume that the two lines 
contribute roughly equally t o the total line strength such 
as observed in GRB 050730 (|Prochaska et al.ll2007af) . 

The majority of the singly ionized species listed in 
Table [5] have ionization potentials above 13.6 eV, and 
with the high column densities these are generally the 
dominant species in DLAs. The metallicity of the dif- 
ferent elements are given relative to the solar value 
[X/H] = log (Nx/Nr) - log (N x /N H ) e in column 5 in Ta- 
ble \5\ Where column densities from different ionization 
levels are derived, the metallicity corresponds to their 
sum. Different elements show relative metallicities in the 
range of 1%— 10% solar. In high column density media, 
the dominant state of Ca is Cam, so the low relative 
metallicity of Ca n does not reflect the real metallicity of 
this element. Also the Mgn lines are saturated, so the 
metallicity is a lower limit. 

To check the derived column densities for consis- 
tency, we also use the apparent o ptical depth (AOD) 
method (jSavage fe Sembachl I1991D . which is valid for 
non-saturated lines. In low resolution data, intrinsically 
saturated lines are smoothed and do not appear satu- 
rated. In practice we find that lines with residual flux 
< 0.7 in the composite spectrum are most likely to be 
saturated, and using these lines will only give a lower 
limit for the column densities. Column 6 in Table [5] lists 
the column densities of the various species, where the in- 
dividual lines were integrated from -600 km s" 1 to +600 
km s . We avoided as best possible the individual sat- 
urated lines as well as lines which are blended with other 
species. Comparing columns 2 and 6, the agreement be- 
tween column densities derived from the COG and the 
AOD methods is good. The exceptions are Fei and Mgi, 
where the AOD method depends strongly on what lines 
are used to derive the column densities. For example 
the Mgi A2852 suggests a much lower column density 
log JV=13.06±0.10 compared to that derived from the 
other two lines (Mgi A1748 and Mgi A1828 which indi- 
cate log7V=14.20±0.10). The large difference in column 
density can be due to the fact that the Mg II A2852 line is 
saturated. Further, these lines also exhibit a large scat- 
ter in the COG plot, and may as explained above be the 
result of a physical difference in the medium probed by 
the GRBs. 

8. DISCUSSION AND FUTURE PROSPECTS 
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Fig. 7. — Curve-of-growth with the best Doppler parameter (and ±lcr uncertainties) fit estimated from ten Fell lines. Other elements, 
are overplotted and all best fit column densities are listed in Table [5] 



Having analyzed a series of parameters, which may af- 
fect the strengths of the metal absorption lines in GRB 
afterglows we arrive at the following results. 

• Variations in excited states are apparent. 

• No variations in the ground states. 

• No correlation with GRB energy release. 

• No correlation with optical afterglow brightness 
one day after the GRB. 

• Dark bursts have stronger lines. 

These results suggest that the medium surrounding the 
GRB at a significant distance is mainly responsible for 
the strengths of the absorption lines, and that the GRB 
itself plays a minor role. Only the excitation of the fine- 
structure lines by the immediate ionizing radiation ap- 
pear to reflect the impact of the GRB explosion. That 
some afterglows are bright and have cither strong or weak 
absorption lines, may just reflect the diverse nature of the 
ISM along different lines of sight through high redshift 
galaxies. 

Characterizing the physical parameters from the com- 
posite is difficult because each absorption line has a dif- 
ferent sample of individual GRB spectra contributing to 
it. As such the composite is not relevant for determin- 
ing the physical conditions present around a single GRB 
because it is a global average of different lines of sight 



at a range of redshifts, where diverse conditions such as 
metallicities, abundance ratios, and densities may pre- 
vail. The results are therefore not directly applicable to 
a particular GRB environment. Biases may also hamper 
the comparison to individual GRB afterglows observed 
in the future, since the individual spectra which domi- 
nate the composite spectrum are relatively bright, and 
are not affected severely by dust obscuration. 

These effects will inevitably influence the column den- 
sities of species, where their individual absorption lines 
lie far apart in wavelength. Although we estimate metal- 
licities from several elements to be 1%— 10% solar, we 
do not attempt to analyze the abundance patterns in 
comparison to the solar or depletion onto dust grains, 
because of the inherent uncertainty of the GRB sam- 
ple. Interestingly, previous studies of individual GRB 
afterglows observed with low resolution spectroscopy 
have indicated ISM metallicities from 10%-100% solar 
dSavaglioll2006t iFvnbo et alj|2006t IVreeswiik et al.ll2006t 



Berger et all 120061; iProchaska et al.ll2007bt iD'Elia et all 
20071 ). althou gh a very low metallicity is reported for 
GRB 090926 (|Rau et alJl20Toh . Conversely, afterglows 
observed at high spectral resolution show relatively lowe r 
metallicities (average 2% solar in iLedoux et al.l [2009( 1. 
The composite low resolution spectrum presented here 
is in better agreement with the high-resolution results 
and the more metal poor systems observed at low reso- 
lution. For multi-phased absorption systems the column 
densities and metallicities can be underestimated using a 
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Fig. 7. — continued. Other elements are fit to the Fen COG. 

COG analysis of low spectral resolution data ()Prochaskal 
2006), which would then cause an even stronger discrep- 
ancy between the high- and low-resolution samples. 

Although the composite spectrum cannot be used 
as an absolute reference scale to measure the relative 
strengths of a set of absorption lines, it can be used 
to discriminate special bursts which truly have outly- 
ing properties such as the pure high-ionization lines i n 
GRB 090426 (ILevesaue et all [MoT IThone et al.1 f20Toh . 
With the recent advent of the wide wavelength coverage 
(3000-25 00 A), second generati on VLT spectrograph, 
X-shooter (|D'Odorico et al.ll2006l ). this composite spec- 
tru m can serve as a reference for li ne identifications (see 
e.g. Jde Ugarte Postigo et al"1l2010[) in future GRB after- 
glow studies. A similar composite spectrum made from 
higher spectral resolution data would help to disentangle 
any of the temporal variations in the strong blends of 
lines from fine-structure and non-excited transitions. 

The composite spectrum can also be used to derive the 
rcdshifts of very noisy spectra using a cross-correlation 
technique. We attempted this for GRB 060614 which did 
not show any detectable absorption lines in its spectrum 
(F09), but fail to determine a reliable redshift probably 
because of the absence of absorption lines rather than a 
bad S/N ratio. On the other hand, we determine tenta- 
tive redshifts for GRB 070129 (z = 2.985), GRB 070508 
(z = 2.790), and GRB 080523 (z = 2.050) consistent 
with possible, albeit noisy detections of either Siiv, Civ 



or Al in lines in their spectra, and also consistent with 
the upper redshift limits in F09. 

An important future study is to investigate the ef- 
fect of the location of the GRB with respect to the host 
galaxy. The locations of high redshift GRBs follow the 
light distributions, and are typically found close to the 
brightest part of the ir host galaxy jBloom et al.l 120021 : 
iFruchter et all 120061: ISvensson etail l2010h . We do not 
have the information of the spatial offsets available to 
make a comparison here. Logically, one could expect 
stronger absorption lines arising from GRBs which are 
located at smaller impact parameters to their hosts, but 
this hypothesis has to be verified first. Another line of 
investigation would be the correlation of absorption line 
strength versus the host galaxy brightness. Again we do 
not have sufficient data available to make this investiga- 
tion. As discussed in Sect. El LBG galaxies show gener- 
ally stronger lines, but are also selected among the intrin- 
sically brighter high redshift galaxies which presumably 
are more massive and metal-rich than the average GRB 
host. Combined with the brightness, it would also be 
an interesting extension to investigate the correlations 
of absorption line strengths with the SFRs and specific 
SFRs of the GRB hosts. 
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TABLE 1 
Weighted Composite 
Afterglow Spectrum 



* ( A ) fx fx uncertainty 



701 


0.003 


0.027 


702 


0.018 


0.027 


703 


0.000 


0.027 


704 


0.000 


0.027 


705 


0.000 


0.027 


706 


0.026 


0.026 


707 


0.000 


0.026 


708 


0.001 


0.025 


709 


0.001 


0.025 


710 


0.000 


0.025 


711 


0.000 


0.025 


712 


0.007 


0.024 


713 


0.030 


0.024 


714 


0.012 


0.024 



Relative flux in the composite 
spectrum and its associated er- 
ror spectrum. This table is avail- 
able in its entirety in a machine- 
readable form in the online jour- 
nal. A portion is shown here for 
guidance regarding its form and 
content. 



TABLE 2 
Absorption Line List 



Line ID A lab (A) c 



<^obs 



W r (A) 



notes 



Hi 


1215. 


67 


Nv 


1238. 


.82 


Nv 


1242. 


.80 


Sn 


1250. 


.58 


Sn 


1253. 


52 


Sn 


1259. 


.52 


Si n 


1260. 


42 


Si n* 


1264. 


,74 


Ci 


1277. 


25 


Ci* 


1277. 


.28 


Oi 


1302. 


.17 


Si n 


1304. 


37 


Oi* 


1304. 


.86 



1239.01 
1243.01 
1249.95 
1252.91 

1259.95 

1264.58 



73 

0.14±0.02 
0.07±0.02 
0.15±0.02 
0.24±0.02 

1.26±0.02 

0.66±0.02 



1277.51 0.09±0.02 



blends 



1302.97 2.29±0.02 a 



Sin* 
Nin 



1309.28 
1317.22 



1309.32 0.27±0.02 
1318.00 0.11±0.02 



Ci 1328.83 
Ci*,** 1329.09-1329.60 



1328.95 0.08±0.01 



a 

blends 



Cn 
Cn* 



1334.53 
1335.71 



1334.99 1.73±0.02 



Cli 


1347. 


.24 


Nin 


1370. 


.13 


Si iv 


1393. 


76 


Si iv 


1402. 


.77 


Gail 


1414. 


40 


Ni n 


1415. 


.72 


CO 


1419. 





Ni n 


1454. 


,84 


Znl 


1457. 


,57 



1347. 


01 





.20±0 


.02 


b 


1370. 


51 





13±0. 


.02 




1393. 


51 





95±0. 


.02 




1402. 


61 


0. 


68±0 


02 




1414. 


95 





16±0 


.02 


a, b 












a, b 


1419. 


00 





06±0 


01 


b 


1456. 


00 





08±0 


02 


a, b 
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Line ID A lab (A) c A obs W r (A) notes 



CO 


1477.5 


1478.4 


0.21±0.02 




Sin 


1526.71 


1526.02 


0.93±0.02 




Sin* 


1533.43 


1533.51 


0.42±0.02 




Co II 


1539.47 


1539.00 


0.03±0.01 


b 


Civ 
Civ 


1548.20 
1550.77 


1548.33 


2.18±0.03 


a 
a 


Ci 


1560.31 


1559.51 


0.09±0.02 


b 


Fen* 


1570.25 


1569.92 


0.08±0.02 


b 


Fc ii* 


1 12 C\ O A C\ 

1602.49 


"i c r\ o cr n 

1602.56 


O OO 1 A Ol 

0.08±0.01 




Fe II 


1608.45 


1 cnc no 
1608.08 


n o c;j.n no 
0.85±0.02 


a 


Fen 


1611.20 


1610.81 


0.18±0.02 


a 


Fen* 


1618.47 


1618.02 


0.03±0.01 


a 


Fen* 


1621.69 


1621.60 


0.11±0.02 


a 


Fe II* 
Fe II* 


1 fion 1 a 
I62y . 16 

1631.13 


1629.50 


0.10±0.02 


a 


Fen* 


1634.35 






a 


Fen* 


1636.33 


1637.55 


0.50±0.02 




Fen* 


1639.40 








Ci 

Ci*,** 


1656.93 

1656.27-1658.12 


1 (\*\7 0.9 


n 1 4-1-0 09 


a 

blcn< 


Alii 


1670.79 


1670.58 


1.04±0.02 




Si i 


1693.29 


1693.03 


0.07±0.02 




Ni ii 


1703.41 


1702.48 


0.14±0.02 




Ni ii 


1709.60 


1709.45 


0.08±0.02 




Nin 


1741.55 


1742.02 


0.14±0.02 




Mgi 


i t a +7 to 

1747.79 


"i t a o no 

1748.02 


n o cr I r\ o 1 

0.05±0.01 


b 


Nin 


1751.92 


"i tct o r\c 

1752. Uo 


r\ oo I r\ oo 




Ni ii 


1804.47 


1805.31 


0.03±0.01 


a 


Si 
Si ii 


1807.31 
1808.01 


1 SOS 1 A 




a 
a 


Si ii* 
Sin* 


1816.93 
1817.45 


1816.95 


0.12±0.01 


a, b 
a, b 


Mgi 


1827.93 


1827.87 


0.09±0.02 




Nin 


1842.89 


1844.27 


0.12±0.02 


b 


Aim 


1854.72 


1854.83 


0.89±0.02 




Aim 


1862.79 


1863.56 


0.68±0.02 




Fel 


1875.16 


1 O 7r A A r* 

1874.46 


r\ "i r\ 1 o o 1 

0.10±0.01 


a 


Fe I 


1 O O O TO 

1883.78 


1 O OO C~7 

1882.67 


A 1 7 I A OO 

0.17±0.02 


a 


Sim 


1892.03 


"1 O OO P" o 

1890.58 


O "I O I o o o 

0.10±0.02 


b 


Fen 


1901.77 


1901.36 


0.14±0.02 


a 


Ti II 
Tin 


1905.77 
1906.24 


1905.99 


0.11±0.02 


a 


Ti II 


1910.75 


1910.56 


0.14±0.02 


a 


Fel 


1937.27 


1938.04 


0.09±0.02 


a 


Co II 


1941.29 


1941.97 


0.07±0.01 


a, b 


TT 

Oo II 


2012.17 


OO "1 o oo 

2010.99 


r\ nc i r\ oo 

0.06±0.02 


b 


Cr II 


2017.57 


2017.61 


0.08±0.02 


b 


Zn II 


2026.14 








Cm 


2026.27 


2025.97 


0.60±0.02 


a 


Mgi 


2026.48 








Cr ii 


2056.26 


2055.51 


0.19±0.02 




Cr ii 
Znll 


2062.23 
2062.66 


2062.67 


0.53±0.02 


a 
a 
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Line 1L) 


Alab (A) 


A obs 




notes 


Or TT 

or 11 


ZUOO. ID 




n 1 9-t-fl 09 
U. 1ZICU.UZ 


a 


Nill* 
Fel 


2166.23 
2167.45 


2167.65 


0.26±0.02 


a 


Nill* 


2175.22 


2175.83 


0.07±0.01 




Mm 


2185.59 


2186.96 


0.51±0.02 




Nin* 


2217.2 


2217.02 


0.24±0.01 




Ni II 


2223.0 


2224.98 


0.07±0.01 




rC 11 


99/1G fiQ 

zz4y.oo 


oo^n no 


n qi _i_n no 




Fen 


2260.78 


2261.00 


0.38±0.02 




Ni II 


2316.7 


2316.36 


0.19±0.02 




rC 11 


OQOQ 1 1 
ZoZo. 1 1 


OQ07 at. 
ZoZ ( .DO 


n 1 0-1-0 09 




Fen* 


2333.52 


2333.30 


0.34±0.02 




ThYi TT 

r u 11 

T?« TT* 

re 11 


OQAA 91 
zo4D.UU 


2346.01 


1.74±0.02 


a 


Fen* 


2349.02 


2349.60 


0.26±0.02 


a 


Fen* 


2359.83 


2360.01 


0.28±0.01 




Fen* 


2365.55 


2365.40 


0.18±0.01 




Fen 


2374.46 


2374.46 


1.00±0.02 




Fen* 


2381.49 






a 


Fen 


2382.77 


2383.59 


1.65±0.02 


a 


Fen* 


2383.79 






a 


Fen* 


2389.36 


2389.40 


0.18±0.02 


a 


PV, TT* 

rC 11 


9QQfi 1 C, 

zoyo. 10 






a 


Fen* 


2396.36 


2396.50 


0.71±0.02 


a 


Fen* 


2399.98 






a 


Fen* 
Fen* 


2405.16 
2407.39 


9406 04 


n 40+0 02 


a 
a 


Fen* 


2411.25 






a 


Fen* 


2411.80 


2411.68 


0.51±0.03 


a 


Fen* 


2414.05 






a 


Sc II 
Sc II 


2561.00 
2563.97 


2562 52 


n i s+o 01 


a 
a 


Mno 


2576.88 


2577.04 


0.45±0.02 




Fen 


2586.65 


2586.49 


1.33±0.02 


a 


Mnll 


2594.50 


2594.36 


0.45±0.02 




Fen* 
Fen 


2599.15 
2600.17 


2600.23 


1.85±0.03 


a 
a 


Mnll 
Fe II* 


2606.46 
2607.87 


2607.53 


0.56±0.01 


a 
a 


P TT* 

re 11 

T?a TT* 

re 11 


9^1 9 fi^ 
ZD1Z.DO 

Z014.D1 


2613.11 


0.51±0.01 


a 


Fen* 


2618.40 


2618.51 


0.06±0.01 


a 


Fen* 


2621.19 






a 


Fell* 


2622.45 






a 


Fen* 
Fell* 


2626.45 
2629.08 


2629.82 


0.90±0.02 


a 
a 


Fe II* 


2631.83 






a 


Fen* 


2632.11 






a 


Fen* 


2740.4 


2739.45 


0.07±0.01 


b 


Fen* 


2747.9 


2749.50 


0.16±0.01 


b 


Fen* 


2756.28 


2756.50 


0.08±0.01 


b 


Mgn 


2796.35 


2796.21 


1.71±0.02 




Mgn 


2803.53 


2803.50 


1.47±0.02 




Mgi 


2852.96 


2852.97 


0.78±0.01 
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Line ID A lab (A) c A obs W r (A) notes 



Fel 


2967.77 


2966.24 


0.05±0.01 


b 


Fel 


2984.44 


2983.49 


0.05±0.01 




Ti II 


3073.88 


3076.01 


0.08±0.01 




Ti II 


3242.93 


3239.56 


O.lliO.Ol 




Ti II 


3384.74 


3385.85 


0.15±0.02 


b 


Can 


3934.78 


3933.97 


0.76±0.02 




Can 


3969.59 


3969.98 


0.66±0.02 




Cal 


4227.92 


4226.93 


0.11±0.02 


b 


MgH 


5209.45 


5210.75 


0.09±0.02 


b 



Lines are blended. 

Line only present in one of the two subsamples, when splitting the full sample by signal-to-noise. 
Vacuum wavelengths. 



TABLE 3 

Absorption line list blueward of Lyck 



Line ID 


Alab (A) 


Aobs (A) 


notes 


Sin 


1206.500 


1206.99 




Ni 
Mnn 


1199.549-1200.710 
1199.391 


1200.47 


blend 


Sin 

Ci*,** 


1190.4158, 1193.2897 
1190.254-1194.686 


1 1 no nfi 


blend 


Ci*,** 


1157.910-1160.51 


1160.52 


a 


Fen 


1142.3656-1144.9379 


1143.00 


a 


Q j* ** 


1138.384-1141.678 


1138.20 


a 


Fen 


1125.4478 


1127.01 




C i*,** 


1121.453 - 1123.460 


1121.49 


a 


Fen 
Nil 


1081.8748 
1083.990 


1085.04 




An 


1066.660 






Fen 


1062.152-1063.9718 


1066.00 




Siv 


1062.662 






An 


1048.2199 


1049.58 




Ovi 
Cii 


1037.6167 
1036.3367 


1037.55 




Ovi 


1031.9261 


1033.00 




Hi 

Mgn 


1025.7225 (Ly/3) 
1025.9681,1026.1134 


1024.39 


a 


Sin 


1020.6989 






Sin 


1012.502 


1013.04 




Nn 
Sin 


989.799 
989.8731 


989.49 




Cm 


977.020 


979.08 




Hi 


972.5368 


973.03 




Ni 


953.415,953.655 


954.49 




Hi 


949.7431 


950.52 




Hi 


937.8035 


937.02 




Hi 


930.7483 


930.99 





Lines arc blended. 
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TABLE 4 

GRB SPECTRA CONTRIBUTING AT LYMAN LIMIT 



GRB name redshift AA r (A) a S/N (A < 912.(1 + z)) logAr(Hi c m - 2 ) b 



GRB 050730 


3.9692 


800- 


-900 





.44 


22 


.10±0. 


10 


GRB 050908 


3.3467 


840- 


-900 


1 


.00 


17 


.60±0. 


10 


GRB 060206 


4.0559 


800- 


-900 





.17 


20 


.85±0. 


10 


GRB 060210 


3.9133 


800- 


-900 





.02 


21 


.55±0. 


15 


GRB 060526 


3.2213 


845- 


-900 





.12 


20 


.00±0. 


15 


GRB 060707 


3.4240 


810- 


-900 





.08 


21 


.00±0. 


20 


GRB 060906 


3.6856 


878- 


-900 





.10 


21 


.85±0. 


10 


GRB 060927 


5.4635 


800- 


-900 





.10 


22 


.50±0. 


15 


GRB 061110B 


3.4344 


873- 


-900 





.52 


22 


.35±0. 


10 


GRB 080603 


2.6892 


870- 


-900 





.12 


21 


.85±0. 


05 


GRB 080607 


3.0368 


800- 


-900 





.04 


22 


.70±0. 


15 



a Wavelength range to which each spectrum contributes to the composite blueward of the Lyman limit, 
k H I column densities are adopted from F09. 



TABLE 5 
Column Densities 



Element 


logA^cm" 2 ) (COG) 


log(X/H) 


log(X/H) a 


[X/H] 


logA^cm" 2 ) (AOD) 


Hi 


22.00±0.10 












Ci 


13.73±0.15 


-8.27±0. 


.18 


-3.57 




13.72±0.09 


Cn 


16.76±0.12 b 


-5.24±0. 


.16 




-1.55 c 


>15.0 


Civ 


16.25±0.15 


-5.75±0. 


.18 






>15.0 


Nv 


13.88±0.05 


-8.12±0. 


.11 


-4.17 




13.90±0.10 


Mgn 


14.91±0.06 


-7.09±0. 


.12 


-4.40 


>-2.68 c 


>14.0 


Mgi 


13.06±0.10 


-8.94±0. 


.14 






14.20±0.15 


Alo 


14.80±0.20 d 


-7.20±0. 


.22 


-5.55 


-1.46 c 


>13.4 


Aim 


14.54±0.20 


-7.46±0. 


.22 






>14.0 


Sin 


16.07±0.22 c 


-5.93±0. 


.24 


-4.49 


-1.27 c 


15.74±0.06 


Si iv 


15.50±0.20 


-6.50±0. 


.22 






>14.1 


Sn 


15.37±0.05 


-6.63±0. 


.11 


-4.88 


-1.75 


15.35±0.12 


Can 


13.17±0.12 f 


-8.83±0. 


.16 


-5.66 


>-3.17 


13.05±0.15 


Ti II 


13.03±0.25 


-8.97±0. 


.27 


-7.05 


-1.92 


12.77±0.23 


Cm 


13.77±0.05 


-8.23±0. 


.11 


-6.36 


-1.87 


13.68±0.09 


Mno 


13.56±0.08 


-8.44±0. 


.13 


-6.57 


-1.87 


13.37±0.03 


Fen 


15.70±0.07 


-6.30±0. 


.12 


-4.50 


-1.80 c 


15.49±0.10 


Fei 


13.50±0.15 g 


-8.50±0. 


.18 






14.90±0.10 


Con 


13.83±0.13 


-8.17±0. 


.16 


-7.01 


-1.16 


13.73±0.21 


Nin 


14.09±0.19 


-7.91±0. 


.21 


-5.78 


-2.13 


14.22±0.21 


Znn 


13.75±0.15 


-8.25±0. 


.18 


-7.44 


-0.81 


13.64±0.04 



a Solar photosphere abundances of the neutral species of each element from Asplund et al. (2009). 
k Blended with C II* A1335, assuming the two lines contribute equally to the measured W r . 
c Derived from the sum of the ionization levels. Lines are saturated, so a lower limit can be determined. 
^ Determined from the single line Alu A1670. 
Contaminated by Sill*. 
The dominant species is Cam. 

s Wide wavelength range for Fei (1875—2984 A). Hence the composite probe different regions, making the fit to the curve of growth bad, 
and values between 13.5 < logA r (Fel) < 14.5 are possible. 



